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We firstly examine hadroproduction of prompt J ftp's at the Fermilab Tevatron in a Monte Carlo framework 
by means of the event generator PYTHIA 5.7 in which those colour-octet processes relevant for charmonium 
production have been implemented accordingly. We find that colour-octet matrix elements presented in literature 
from pp collider data are systematically overestimated due to overlooking of the effective primordial transverse 
momentum of partons (i.e. including higher-order QCD effects). We estimate the size of these effects using 
different parton distribution functions. Finally, after normalization to Tevatron data, we extrapolate up to LHC 
energies making a prediction on the expected pt differential cross-section for charmonium. 



The experimental discovery at Fermilab jy of 
an excess of inclusive production of prompt heavy 
quarkonia (mainly for J/ip and ijj(2S) resonances) 
in antiproton-proton collisions triggered an in- 
tense theoretical activity beyond what was con- 
sidered conventional wisdom just few years ago 
|2| . Indeed, the discrepancy between the so-called 
colour-singlet model (CSM) in hadroproduction 
|3| and the experimental data amounts to more 
than an order of magnitude and cannot be at- 
tributed to those theoretical uncertainties arising 
from the ambiguities on the choice of a particu- 
lar parton distribution function (PDF), the heavy 
quark mass or different energy scales. 

Recently it has been argued that the surplus 
of charmonia production can be accounted for by 
assuming that the heavy quark pair not necessar- 
ily has to be produced in a colour-singlet state at 
the short-distance partonic process itself [f| . Con- 
versely, it can be produced in a colour-octet state 
evolving non-perturbatively into quarkonium in a 
specific final state with the correct quantum num- 
bers according to some computable probabilities 
governed by the internal velocity of the heavy 
quark. This mechanism, usually named as the 
colour-octet model (COM), can be cast into the 
rigorous framework of an effective non-relativistic 
theory for the strong interactions (NRQCD) de- 
riving from first principles ||. 

However, the weakness of the COM lies in the 
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fact that the non-perturbative parameters charac- 
terizing the long-distance hadronization process 
beyond the colour-singlet contribution (i.e. the 
colour-octet matrix elements) are so far almost 
free parameters to be adjusted from the fit to ex- 
perimental data, though expected to be mutually 
consistent according to the NRQCD power count- 
ing rules. 

On the other hand, an attractive feature of the 
colour-octet hypothesis consists of the universal- 
ity of the NRQCD matrix elements entering in 
other charmonium production processes. Let us 
look below in some detail at the way hadroniza- 
tion is folded with the partonic description of 
hadrons, focusing for concreteness on the couple 
of related papers || \7\. 

In the first paper [3|, Cho and Leibovich con- 
sidered for charmonium production only the 3 <Si 
coloured intermediate state as a first approach, 
computing the matrix elements as products of 
perturbative parts for the short-distance partonic 
processes, and the colour-octet matrix element 
concerning the long-distance hadronization, in a 
similar but generalized way as in the CSM. Fi- 
nally, a convolution of concrete parton distribu- 
tion functions and the differential cross-section 
for the QQ production subprocess was performed, 
whereby the pt dependence of the charmonia pro- 
duction coming from the latter. Moreover, in ref. 
the same authors take into account further 
contribution from new coloured states (for more 
details see the quoted references) concluding fi- 
nally that at high enough p t a two-parameter fit 
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is actually required to explain the observed inclu- 
sive pt distribution of charmonia production at 
the Tevatron. 

However, it is well-known for a long time that 
higher-order effects (K factors) play an important 
role in inclusive hadroproduction ||. Namely, 
beyond the primordial transverse momentum kt 
of partons in hadrons, higher-order QCD effects 
such as initial-state radiation of gluons by the in- 
teracting partons add up to yield an effective in- 
trinsic transverse momentum which certainly has 
to be considered in charm production even at high 
Pt [||. As we shall see, if overlooked at all, the ef- 
fect on the fit parameters (and ultimately on the 
colour-octet matrix elements) amounts to a sys- 
tematic overestimate. 




Figure 1. Curves obtained from PYTH1A (not 
fit) including colour-octet mechanism for prompt 
J/ip production at the Tevatron using the same 
parameters as in Ref. [7]. The solid line corre- 
sponds to initial- and final-state radiation turned 
off and the dotted line to initial- and final-state 
radiation on. The MRSDO parton distribution 
function was employed as in [7]. 



In this work we have implemented in the event 
generator PYTH1A 5.7 [|o| those new colour- 
octet mechanisms for the dominant gluon-gluon 



fusion process 0. Let us stress that in our analy- 
sis initial-state radiation was incorporated within 
the framework of the Lund Monte Carlo frame- 
work (If) 0. 
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Figure 2. Two-parameter fit to the experimental 
Tevatron data, where initial- and final-state radi- 
ation were incorporated via PYTH1A generation. 
The CTEQ 2L parton distribution function was 
employed in this plot. 



Using the same numerical values for the colour- 
octet matrix elements as reported in table I of 
Ref. if initial- and final-state radiation are 
turned off, a relative good agreement of the the- 
oretical curve and the experimental points is ob- 
tained (see Fig. f), as should be reasonably ex- 
pected 0. However, if the initial- and final-state 



^Originally PYTHIA generates direct production of J ftp's 
via the CSM Jip| . Technically, we have implemented in the 
PYTHIA subroutine PYSIGH the two extra differential 
cross-sections for the process g + g^J/ip + g corresponding 
to the coloured 3 Si and 1 Sn + 3 Po states whose matrix 
elements can be found in flfJ M . 

4 Initial-state radiation is switched on bv default in 
PYTHIA: MSTP(61) is themaster switch 0]. 
5 The same PDF as in Ref. jjj (MRSDO) was used. Still the 
theoretical curve stands slightly below the experimental 
points maybe due to the fact that we have neglected other 
non-dominant channels like qq^J/ipg and qg^J/ipg |7| 
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Table 1 

Colour-octet matrix elements (in units of GeV 3 ) from the fit of Tevatron data on prompt J/ip produc- 
tion for different parton distribution functions. For comparison we quote the values given in Ref. [7]: 
(6.6±2.1)xl(T 3 and (2.2±0.5) x 1(T 2 respectively. 



matrix element 
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Figure 3. Our prediction for prompt J/ip direct 
production at the LHC using PYTHIA with the 
new colour-octet matrix elements a) dotted line: 
CTEQ 2L, b) solid line: GRV HO, and c) dot- 
dashed line: MRSDO. The rapidity cut \y\ < 2.5 
on the J ftp was required. 



radiation ^ are switched on the predicted curve 
stands well above the experimental data. By in- 
spection we observe from Fig. 1 that the the- 
oretical curve including initial-state radiation is 
shifted upwards with respect to the experimen- 

fact that we somehow underestimate the production rate 
of prompt charmonia reinforces our arguments favouring 
lower values for the colour-octet matrix elements. 
6 It should be noted that initial-state radiation and final- 
state radiation have opposite effects in the pt spectrum, 
the former enhancing the high pt tail whereas the latter 
smears the distribution. 



tal distribution, in accordance with the expected 
"kick" caused by the effective primordial trans- 
verse momentum of partons || ||. Accordingly, 
keeping radiation effects it turns out that the val- 
ues for the colour-octet matrix elements have to 
be lowered by significant factors (!). Of course, 
the PYTHIA treatment of the effective kt is not 
guaranteed to be perfect but, nevertheless, should 
give a reasonable estimate of such effects. 

In order to assess the importance of the effec- 
tive intrinsic transverse momentum of partons we 
have made three different choices for the proton 
PDF []: 

a) the leading order CTEQ 2L (by default in 

PYTHIA 5.7) 

b) the next to leading order MRSDO (the same 

as used in j?J) 

c) the next to leading order GRV 94 HO 

As already mentioned before, the theoretical 
curve of the inclusive p t distribution of prompt 
J/^'s stands in all cases above Tevatron exper- 
imental points if the set of parameters from J?]] 
are blindly employed in the PYTHIA generation 
with initial-state radiation on. Motivated by this 
systematic discrepancy, we performed new fits 
for the prompt J/tp direct production at Teva- 
tron (feed-down from radiative decay of \cj reso- 
nances was experimentally removed). The corre- 
sponding colour-octet matrix elements are shown 
in table I and the plot for case a) in Fig. 2. 

7 See Jf2|] for technical details about the package of Par- 
ton Density Functions available at the CERN Program 
Library. References therein. 
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Finally, we have generated prompt J ftp's in 
proton-proton collisions at LHC energies (center- 
of-mass energy = 14 TeV) by means of our 
"modified" version of PYTHIA with the colour- 
octet matrix elements as shown in Table I - i.e. 
after normalization to Tevatron data. Indeed, an 
"order of magnitude" estimate of the expected 
production of charmonia at the LHC is suitable 
from many points of views Jl3| : J/ip can be con- 
sidered either as a signal in its own right or as 
a source of background for other interesting pro- 
cesses involving J/tp's, like CP studies from the 
cascade channel — > J/tp K®. In Fig. 3 we show 
the pt inclusive distributions for direct prompt 
J/ip's at the LHC obtained for each of the three 
PDF's employed in our study. Such predictions 
should be taken, however, with a grain of salt 
since gluons emitted either at the initial or at the 
final state could contribute at high energies to 
the final yield of charmonia via the colour-octet 
mechanism. (This possibility is currently being 
investigated.) 

In summary, we have investigated higher-order 
effects induced by initial-state radiation on the 
extraction of the NRQCD matrix elements from 
hadroproduction at high p t by means of an event 
generator (PYTHIA 5.7) with colour-octet mech- 
anisms implemented in. We conclude that for 
different PDF's the overlooking of the effective 
primordial kt leads sytematically to a significant 
overestimate of the colour parameters. We have 
derived new colour-octet matrix elements for J/ ip 
production from Tevatron data using three differ- 
ent sets of PDF's. 

Our conclusion is that a more detailed analy- 
sis is actually required to extract NRQCD colour- 
octet matrix elements from hadroproduction data 
in a reliable way thereby allowing a meaning- 
ful comparison with other charmonia production 
processes, like photoproduction at HERA |14) . 
Finally we have estimated the prompt J/tp in- 
clusive direct production at the LHC for three 
different PDF's, finding an overall relative good 
accordance at moderate transverse momentum. 
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